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DISPOSITIF DE SURVEILLANCE DE L'INTEGRITE DES INFORMATIONS 
DELIVREES PAR UN SYSTEME HYBRIDE INS/GNSS 

L'invention est relative a la surveillance de I'integrite des 
5 informations de position et de Vitesse issues d'une hybridation entre une 
centrale inertielle et un recepteur de positionnement par satellites. Elle 
concerne plus precisement un equipement de navigation connu dans la 
technique sous le nom de systeme INS/GNSS (de I'anglo-saxon "Inertial 
Navigation System" et "Global Navigation Satellite System") hybride en 
10 boucle fermee, l'hybridation etant dite lache car realisee en axes 
geographiques. 

Une centrale inertielle est constitute d'un ensemble de capteurs 
inertiels (capteurs gyrometriques et capteurs accelerometriques) associe a 
une electronique de traitement. Une plate-forme de calcul, appelee plate- 
15 forme virtuelle PFV, delivre alors les informations de Vitesse et de position du 
porteur dans un referentiel precis (souvent note TGL, Triedre Geographique 
Local). La plate-forme virtuelle PFV permet la projection et ['integration des 
donnees issues des capteurs inertiels. La centrale inertielle fournit des 
informations precises a court terme mais derivant sur le long terme (sous 
20 influence des defauts capteurs). La maTtrise des defauts capteurs 
represente une proportion tres importante du cout de la centrale inertielle. 

Un recepteur de positionnement par satellites fournit des 
informations de position et de vitesse du porteur par triangulation a partir des 
positions de satellites defilants visibles du porteur. Les informations fournies 
25 peuvent etre momentanement indisponibles car le recepteur doit avoir en vue 
directe un minimum de quatre satellites du systeme de positionnement pour 
pouvoirfaire un point. Elles sont en outre d'une precision variable, dependant 
de la geometrie de la constellation a la base de la triangulation, et bruitees 
car reposant sur la reception de signaux de tres faibles niveaux provenant de 
30 satellites eloignes ayant une faible puissance d'emission. Mais elles ne 
souffrent pas de derive a long terme, les positions des satellites defilants sur 
leurs orbites etant connues avec precision sur le long terme. Les bruits et les 
erreurs peuvent etre lies aux systemes satellitaires, au recepteur ou a la 
propagation du signal entre I'emetteur satellitaire et le recepteur de signaux 
35 GNSS. En outre, les donnees satellites peuvent etre erronees par suite de 
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pannes affectant les satellites. Ces donnees non integres doivent alors etre 
reperees pour ne pas fausser la position issue du recepteur GNSS. 

Pour prevenir les pannes satellites et assurer I'integrite des 
mesures GNSS, il est connu d'equiper un recepteur de position nement par 
satellites d'un systeme d'estimation de precision et de disponibilite dit RAIM 
(de I'anglo-saxon "Receiver Autonomous Integrity Monitoring") qui se base 
sur la geometrie de la constellation de satellites utilisee lors de la 
tnangulation et sur 1'evolution previsible a court terme de cette geometrie 
deduite de la connaissance des trajectoires des satellites. Cependant ce 
systeme, purement lie au systeme de localisation par satellites, n'est pas 
appl.cable a la surveillance d'un point de localisation issu d'un systeme 
hybride (INS/GNSS) et ne peut detecter que certains types de pannes en un 
temps donne. 

(-'hybridation consiste a combiner mathematiquement les 
15 informations de position et de vitesse fournies par la centrale inertielle et le 
recepteur de positionnement par satellites pour obtenir des informations de 
position et de vitesse en tirant avantage des deux systemes. Ainsi, la 
precision sur le positionnement fourni par le systeme GNSS permet' de 
maitriser la derive inertielle et les mesures inertielles peu bruitees permettent 
20 de filtrer le bruit sur les mesures du recepteur GNSS. Cette combinaison fait 
tres souvent appel a la technique de filtrage de Kalman. 

Le filtrage de Kalman s'appuie sur les possibles de modelisation 
de revolution de I'etat d'un systeme physique considere dans son 
environnement, au moyen d'une equation dite "d'evolution" (estimation a 
25 priori), et de modelisation de la relation de dependance existant entre les 
etats du systeme physique considere et les mesures par lesquelles il est 
percu de I'exterieur, au moyen d'une equation dite "d'observation" pour 
permettre un recalage des etats du filtre (estimation a posteriori). Dans un 
filtre de Kalman, la mesure effective ou "vecteur de mesure" permet de 
30 realiser une estimee a posteriori de I'etat du systeme qui soit optimale dans 
le sens ou elle minimise la covariance de I'erreur faite sur cette estimation. 
La partie estimateur du filtre genere des estimees a priori du vecteur d'etat 
du systeme en utilisant I'ecart constate entre le vecteur de mesure effectif et 
sa prediction a priori pour engendrer un terme correctif, appele innovation. 
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Cette innovation est appiiquee a I'estimee a priori du vecteur d'etat du 
systeme et conduit a I'obtention de i'estimee optimale a posteriori. 

Dans le cas d'un systeme INS/GNSS hybride, le filtre de Kalman 
recoit les points de position et de Vitesse fournis par la centrale inertielle et le 
5 recepteur de positionnement par satellites, modelise revolution des erreurs 
de la centrale inertielle et delivre I'estimee a posteriori de ces erreurs qui sert 
a corriger le point de positionnement et de vitesse de la centrale inertielle. 

L'estimation des erreurs de position et de vitesse dues aux 
defauts des capteurs inertiels apparaissant en sortie de la plate-forme 
10 virtuelle PFV de la centrale inertielle est realisee par le filtre de Kalman. La 
correction des erreurs par le biais de leur estimation faite par le filtre de 
Kalman peut alors se faire en entree de la plate-forme virtuelle PFV 
(architecture en boucle fermee) ou en sortie (architecture en boucle ouverte). 

Lorsque les defauts des capteurs de la centrale inertielle ne sont 
15 pas trop importants, il n'est pas necessaire d'appliquer les corrections en 
entree de la plate-forme virtuelle PFV. La moderation du systeme 
(linearisation des equations regissant revolution du systeme), au sein du 
filtre reste valide. L'estimee a posteriori des erreurs de la centrale inertielle 
calculee dans le filtre de Kalman n'est utilisee que pour l'elaboration de 
20 I'estimee optimale de la position et de la vitesse du porteur compte tenu des 
informations de position et de vitesse foumies par la centrale inertielle et par 
le recepteur GNSS. L'hybridation est alors dite en boucle ouverte. 

Lorsque les defauts inertiels sont trop importants la linearisation 
des equations regissant revolution du modele inertiel integre au sein du filtre 
25 de Kalman n'est plus valide. II est done imperatif d'appliquer les corrections a 
la plate-forme virtuelle PFV pour rester dans le domaine lineaire. L'estimee a 
posteriori des erreurs la centrale inertielle calculee dans le filtre de Kalman 
sert non seulement a l'elaboration de l'estimee optimale de la position et de 
la vitesse du porteur mais aussi au recaiage de la centrale inertielle au sein 
30 de la plate-forme virtuelle PFV. L'hybridation est alors dite en boucle fermee. 

L'hybridation peut egalement se faire en observant des 
informations GNSS de natures differentes. Soit on considere les position et 
vitesse du porteur resolues .par le recepteur GNSS. On parle alors 
d'hybridation lache. Soit on considere les informations extraites en amont par 
35 le recepteur GNSS que sont les pseudo-distances et les pseudo-vitesses 
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(grandeurs directement isssues de la mesure du temps de propagation et de 
I'effet Doppler des signaux emis par les satellites en direction du recepteur). 
On parle alors d'hybridation serree. 

Avec un systeme INS/GNSS en boucle fermee ou le point resolu 
5 par le recepteur GNSS est utilise pour recaler les informations provenant de 
la centrale inertielle, il est necessaire de prefer une attention particuliere aux 
defauts affectant les informations fournies par les satellites car le recepteur 
qui les recoit propagera ces defauts a la centrale inertielle en entraTnant un 
mauvais recalage de cette derniere. Le probleme se pose de facon 
io particulierement critique pour assurer I'integrite d'un point hybride INS/GPS. 

Une maniere connue de proceder a la surveillance de I'integrite 
d'un systeme hybride INS/GNSS en boucle fermee est decrite dans le brevet 
americain US 5,583,774. Elle consiste a espacer les recalages d'un temps 
suffisamment long (par exemple 30 minutes) pour qu'un detecteur a base de 
15 filtres de Kalman surveillant revolution des ecarts de mesures de pseudo- 
distance et de pseudo-vitesse par rapport au porteur de chaque satellite 
visible ait pu isoler les satellites defaillants. 

Un autre precede connu pour surveiller I'integrite d'un systeme 
hybride INS/GNSS est decrit dans le brevet americain US 5,923,286. II 
20 consiste a utiliser un equipement RAIM pour autoriser ou non I'hybridation. 
Lorsque I'equipement RAIM signale une perte d'integrite, I'hybridation est 
gelee et le point de position et de Vitesse est foumi par la centrale INS en 
tenant compte de ses derives et biais mesurees juste avant la perte 
d'integrite. Pour que cela marche, il faut que la centrale inertielle n'ait pas ete 
25 polluee par I'erreur de point commise par le recepteur GNSS, ce qui interdit 
de la recaler avec le recepteur GNSS. Le precede est done reserve 
uniquement au systeme hybrides INS/GNSS en boucle ouverte. 

La presente invention a pour but une surveillance de I'integrite des 
informations de position et de Vitesse d'un systeme hybride constitue d'une 
30 centrale inertielle recalee a I'aide d'un recepteur GNSS en utilisant la position 
et la vitesse resolues du porteur issues du recepteur GNSS. 

Elle a pour objet un dispositif de surveillance de I'integrite d'un 
systeme hybride constitue d'une centrale inertielle INS, d'un recepteur de 
positionnement par satellites GNSS operant a partir d'une constellation de N 
35 satellites visibles et d'un filtre de Kalman d'hybridation ayant un vecteur 
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d'etats correspondants aux erreurs du systeme hybride, en particulier, les 
erreurs de biais et de derive residuelles de la centrale inertielle INS, 
observant les ecarts entre les points de positionnement et de Vitesse fournis, 
en coordonnees geographiques, par la centrale inertielle INS et par le 
5 recepteur GNSS, comportant une matrice devolution F modelisant revolution 
des erreurs du systeme hybride, une matrice d'observation H modelisant les 
relations liant le vecteur d'etats et les ecarts observes entre les positions et 
les vitesses delivrees par la centrale inertielle INS et le recepteur GNSS, et 
un gain K minimisant la covariance de I'erreur faite sur I'estimation a 
10 posteriori des erreurs de position et de Vitesse issues de la centrale 
inertielle, et delivrant une estimee a posteriori des erreurs du systeme 
hybride qui sert a recaler la centrale inertielle. Ce systeme de surveillance 
est remarquable en ce que le recepteur GNSS delivre, en plus d'un point de 
position effectue a partir des N satellites de la constellation qu'il a en vue, N- 
15 1 points de position effectues a partir de la constellation des N satellites 
visibles privee, a chaque fois, d'un satellite different et en ce qu'il comporte 
un circuit detecteur de probleme satellite comportant un banc de N filtres 
predicteurs-estimateurs de I'erreur induite par le satellite qui a ete ote lors de 
la resolution d'un point a N-1 satellites, ayant le gain K et la matrice 
20 devolution F du filtre de kalman d'hybridation, observant chacun I'ecart entre 
le point de position, en coordonnees geographiques, delivre par le recepteur 
GNSS en observant les N satellites visibles et I'un des points de position, 
egalement en coordonnees geographiques, delivres par le recepteur GNSS 
en observant N-1 satellites visibles et des circuits de test comparant les etats 
25 des N filtres predicteurs-estimateurs a leurs variances et detectant une 
panne satellite lorsque le test est positif, I'ecart trouve etant superieur a un 
seuil de detection. 

Avantageusement, les circuits de test ont des seuils de detection 
30 de panne satellite qui resultent de tests statistiques tenant compte de la 
covariance associee au type d'erreur de positionnement considere. 

Avantageusement, les circuits de test ont des seuils de detection 
de panne satellite qui resultent de tests statistiques tenant compte de la 



1 er depot 



6 




covariance associee au type d'erreur de positionnement considere et qui 
sont fonction du taux de fausse alarme acceptable pour le test. 

Avantageusement, le systeme de surveillance comporte un circuit 
5 d'inhibition des recalages intercale entre la sortie du filtre de Kalman 
d'hybridation et une entree de recalage de la centrale inertielle, et active par 
un test positif de detection de probleme satellite. 

Avantageusement, les filtres predicteurs-estimateurs du circuit 
10 detecteur de probleme satellite (premier banc de filtre) ont tous la meme 
matrice d'observation. 

Avantageusement, le recepteur GNSS delivre en plus des points 
de position, en coordonnees geographiques, effectues a partir des N 

15 satellites visibles et ceux effectues a partir de N-1 satellites deduits de la 
constellation formee par les N satellites visibles en retirant a chaque fois un 
satellite visible different, Nx(N-1)/2 points de position avec N-2 satellites 
deduits des N satellites visibles en retirant a chaque fois deux satellites 
visibles differents, et le systeme de surveillance comporte un circuit 

20 identificateur de satellite defaillant compose d'un deuxieme banc de Nx(N- 
1)/2 filtres de predicteurs-estimateurs d'ecarts de positionnement ayant le 
gain K et la matrice devolution F du filtre de Kalman d'hybridation, observant 
chacun I'ecart entre un point de position, en coordonnees geographiques, 
delivre par le recepteur GNSS a partir de N-1 satellites parmi les N visibles et 

25 I'un des points de positionnement delivres par le recepteur de 
positionnement par satellites a partir de N-2 satellites parmi les N visibles 
deduite de la constellation specifique de (N-1) satellites visibles privee de I'un 
de ses satellites visibles, les filtres estimateurs-detecteurs pouvant etre 
regroupes par families de N-2 elements en fonction de la constellation 

30 specifique de N-1 satellites visibles prise en compte, des circuits de test 
comparant les Stats des Nx(N-1)/2 filtres predicteurs-estimateurs par rapport 
a leurs variances afin de detecter une eventuelle anomalie et un circuit de 
traitement identifiant, en cas d'anomalie detectee, un satellite fautif comme 
etant le satellite exclu d'un point a n-1 satellites dont la famille de filtres 
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predicteurs-estimateurs est la seule a n'avoir aucun de ses elements 
detectant une anomalie 

Avantageusement, ies filtres predicteurs-estimateurs du circuit 
5 identificateur de satellite defaillant ont tous la meme matrice d'observation. 

Avantageusement, les filtres predicteurs-estimateurs du detecteur 
de probleme satellite et du circuit identificateur de satellite defaillant ont tous 
la meme matrice d'observation. 

10 

D'autres caracteristiques et avantages de ['invention ressortiront 
de la description ci-apres d'un mode de realisation donne a titre d'exemple. 
Cette description sera faite en regard du dessin dans lequel : 
15 - des figures 1 et 2 represented le schema de principe d'un 

systeme hybride de position nement INS/GNSS de type lache 
en boucle fermee, la figure 1 defaillant la centrale inertielle et 
ia figure 2 le filtre de Kaiman, 
- une figure 3 represente un systeme hybride de positionnement 
20 INS/GNSS pourvu d'un circuit detecteur de probleme satellite 

conforme a I'invention, 
une figure 4 represente un systeme hybride de positionnement 
INS/GNSS pourvu d'un circuit detecteur de probleme satellites 
et d'un circuit d'identificateur de satellite fautif conforme a 
25 ('invention, et 

~ une figure 5 detaille la constitution du circuit identificateur de 
. satellite fautif de la figure 4. 

La figure 1 montre ('architecture d'un systeme hybride de 
30 positionnement equipe d'une centrale inertielle 1 et d'un recepteur de 
positionnement par satellites GNSS 2, diten boucle fermee et de type lache 
car il se serf des informations de position (latitude, longitude et altitude) et de 
vitesse (vitesse nord, vitesse est et Vitesse verticale) resolues du porteur 
delivrees par le recepteur GNSS 2 pour le recalage de la centrale inertielle 
35 INS1. 
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Comme montre, la centrale inertielle comporte un ensemble 10 de 
capteurs inertieis (accelerometres et gyrometres) et de circuits de calcul, 
appeles plate forme virtuelle PFV, pour effectuer les calculs de position, 
attitude et vitesses sur les sorties des capteurs inertieis. La plate-forme 

5 virtuelle PFD comporte principalement un ensemble d'integrateurs 11 
operant sur tes signaux Q des gyrometres afin d'en deduire ('attitude du 
porteur, un changeurde repere 12 permettant, a partir de la connaissance de 
rorientation du porteur (attitude) du porteur, d'exprimer dans un repere 
geographique lie a ia terre, les accelerations acc mesurees par les 

10 accelerometres dans un repere lie au porteur et deux ensembles integrateurs 
successifs 13, 14 operant sur les composantes de ('acceleration selon le 
repere geographique pour en deduire les composantes de la vitesse du 
porteur par rapport au nord Vn, par rapport a Test Ve et par rapport a la 
verticale Vv ainsi que la position du porteur en latitude, longitude et altitude. 

15 Les capteurs inertieis souffrent de defauts inherents a leur conception : des 
biais et des derives qui font que les informations de position et d'attitude 
delivrees par une centrale inertielle se degradent avec le temps de 
fonctionnement La lutte centre ces defauts est couteuse, ainsi il existe 
differentes Glasses de centrale inertielle en fonction des tolerances de biais 

20 et de derives acceptees. 

Le recepteur GNSS 2 est, par exemple un recepteur GPS. II 
mesure les pseudo-distances et pseudo-vitesses du porteur par rapport a au 
moins quatre satellites visibles du porteur du systeme hybride de 
positionnement et resout, par triangulation, la position du porteur en latitude, 

25 longitude et altitude ainsi que les composantes de sa vitesse par rapport au 
nord Vn, par rapport a Pest Ve et par rapport a la verticale Vv. II ne sera pas 
detaille plus avant car il est bien connu de Fhomme du metier et largement 
decrit dans la litterature ou de nombreux ouvrages lui ont ete consacres, par 
exemple le livre de Elliot D. KAPLAN intitule "Understanding GPS Principles 

30 and Applications" paru chez Artech House Publishers ISBN 0-89006-793-7. 
Les informations de position et de vitesse delivrees par un recepteur GNSS 
ne souffrent pas de derive puisqu'eiles proviennent de mesures instantanees 
des eioignements et vitesses d'eloignement de satellites dont les orbites et 
les positions sur leurs orbites sont connues a long terme. Par centre, elles 

35 souffrent d'indisponibilites, lorsque qu'il n'y a pas un minimum de quatre 
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satellites visibles. Elles sont bruitees car elles proviennent du traitement de 
signaux recus qui sont de tres faibles niveaux d'energie en raison de la 
distance d'eloignement des satellites et de leurs faibles puissances 
d'emission et qui peuvent avoir ete corrompus en cours de transmission 

5 (rencontre de particules chargees entre le satellite et le recepteur). Leur 
fiabilite depend de celle des informations transmises par les satellites. 

Par rapport au recepteur GNSS, la centrale inertielle a I'avantage 
de ne pas faire appel a la collaboration d'equipements exterieurs au porteur 
et done de ne pas etre aussi sensible a I'exterieur. L'integrite et la 

10 disponibilite des informations inertielles etant bien meilleures, la centrale 
inertielle conserve tout son interet. Pour lutter contre sa derive et ses biais 
autrement que par des mesures couteuses prises au niveau de ses capteurs 
inertiels, il a ete propose de la recaler periodiquement au moyen du point de 
position et de vitesse fourni par un recepteur GNSS, apres un filtrage 

1 5 antibruit. 

Ce recalage se fait au moyen d'un filtre de Kalman 3 dit 
d'hybridation montre de maniere plus detaillee a la figure 2. 

Un filtre de Kalman permet d'obtenir une estimation de I'etat a 
posteriori d'un systeme, optimale dans le sens ou la covariance de I'erreur 
20 faite sur cette estimation a posteriori est minimale (on parle aussi 
d'estimation optimale au sens des moindres carres). II se base sur une 
modelisation de revolution du systeme et sur une modelisation de la relation 
existant entre la mesure et I'etat du systeme. 

Dans sa plus grande generalite, la modelisation de 1'evolution de 
25 I'etat du systeme est definie par une equation differentielle vectorielle 
linearisee au premier ordre qui s'exprime, apres discretisation du modele 
continu, par une equation devolution de la forme: 

= F K X K + L K U i + W K 

x etant le vecteur d'etat de dimension p, u un vecteur de controle, F K la 
30 matrice d'evolution definissant la relation entre le vecteur d'etat a I'etape K et 
le vecteur d'etat a I'etape K+1 en I'absence de vecteur de controle et de bruit 
de fonctionnement, L K une matrice definissant la relation entre le vecteur de 
contr6le et le vecteur d'etat au cours d'une meme etape et w K un bruit de 
fonctionnement au cours d'une etape, suppose blanc et gaussien a valeur 
35 moyenne nulle. 
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La modelisation de la relation existant entre la mesure et Tetat du 
systeme est definie par une autre equation differentielie du premier ordre qui 
s'exprime, apres discretisation du modele continu, par une equation 
d'observation de la forme: 

5 z K =H K x K +v K 

z etant le vecteur de mesure (grandeurs du systemes physiquement 
observables) de dimension m, Hk la matrice d'observation definissant (a 
relation entre le vecteur de mesure et le vecteur d'etat au cours d'une meme 
etape et vk un bruit de mesure au cours d'une etape suppose bianc et 

10 gaussien a valeur moyenne nulle. 

Le filtre de Kalman montre a la figure 2 ne comporte pas de 
vecteur de controle et done pas de matrice Lk. II est recursif et repose sur la 
determination d'une estimation a priori x KiK _ x du vecteur d'etat du systeme a 
I'etape K a partir de la connaissance a posteriori du vecteur d'etat du 

15 systeme x K _ UK ^ a I'etape precedente K-1 et de ('application a I'estimee a 
priori x K/K _ x d'un terme correcteur dependant de l'ecart constate entre le 
vecteur mesure z K constate au cours de cette etape K et son estimee a 
priori z KfK ^ deduite de 1'estimation a priori x KfK _ x du vecteur d'etat. 

L'estimee a priori z K/K _ x du vecteur de mesure a I'etape K est 

20 determinee par application des equations d'evolution et d'observation a 
1'estimee a posteriorU^,^ du vecteur d'etat correspondant a I'etape K-1 
precedente. Cette operation est illustree sur la figure 2 par la presence d'une 
boucle de retroaction constitute d'un circuit a retard 30 et de deux fiitres 31, 
32. 

25 Le circuit a retard 30 connecte en sortie du filtre de Kalman 

permet de memoriser Testimee a posteriori x K _ VK _ x du vecteur d'etat 
disponible a la sortie du filtre de Kalman au cours de i'etape precedente K-1. 

Le filtre 31 permet, par mise en ceuvre de Pequation d'evolution, 
I'obtention de I'estimee a priori x K/K _ x 6u vecteur d'etat a I'etape K, a partir de 

30 i'estimee a posteriori x K _ UK „ x du vecteur d'etat du systeme a I'etape 
precedente K-1 . Sa fonction de transfert est definie par la matrice de 
I'equation d'evolution. 

Cette estimee a priori x K/K ^ du vecteur d'etat a i'etape K est 
ensuite utilisee pour obtenir, au moyen d'un deuxieme filtre 32 a I'estimee a 

35 priori z K/K ^ du vecteur de mesure a I'etape K par application de I'equation 
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d'observation. Ce deuxieme filtre 32 a, pour ce faire, une fonction de transfert 
definie par la matrice H K de I'equation d'observation. 

L'estimee a priori i KIK _ x du vecteur de mesure a I'etape K obtenue 
dans la boucle de retroaction est appliquee, en entree du fiitre de Kalman, a 

5 un circuit soustracteur 33 qui recoit par ailleurs le vecteur de mesure z K 
effectivement mesure au cours de I'etape K et qui delivre un vecteur d'erreur 
r K , egalement appele innovation, attestant de I'erreur commise lors de 
i'estimation a priori. Ce vecteur d'erreur r K est transforme par un troisieme 
filtre 34 en un vecteur de correction. Ce vecteur de correction est additionne 

10 par un deuxieme sommateur 35 a l'estimee a priori x K/K _ x du vecteur d'etat 
pour I'etape K issue du premier filtre 31, pour obtenir l'estimee a posteriori 
x K/K du vecteur d'etat qui constitue ia sortie du filtre de Kalman. 

Le troisieme filtre 34, qui foumit le terme correctif, est connu sous 
le nom de filtre de gain de recalage. II a une fonction de transfert definie par 

15 une matrice K K determinee de facon a minimiser la covariance de I'erreur 
faite sur I'estimation a posteriori. 

Kalman a montre que la matrice de gain optimal K K pouvait etre 
determinee par une methode recursive a partir des equations: 

- de la matrice de covariance de l'estimee a priori du vecteur 
20 d'etat 

- de definition de la matrice de gain elle-meme 

K K = P K i K -iH k {h k P K ik-\H T k + & K ) 
-de mise a jour de la matrice de covariance de l'estimee a 
25 posteriori du vecteur d'etat 

p K1K = {l~ KkHr ) P KIK-\ 

P etant la matrice de covariance du vecteur d'etat, soit estimee a priori pour 
I'etape K a partir de I'etape K-1 si P est affecte de I'indice K/K-1, soit estimee 
a posteriori si P est affecte de I'indice K-1 /K-1, 

30 R etant la matrice de covariance des bruits d'observation w K , 

Q etant la matrice de covariance des bruits d'etat (ou de fonctionnement) v K . 

A ['initialisation, la matrice de covariance du vecteur d'etat et le 
vecteur d'etat sont pris egaux a leurs estimees les plus vraisemblables. Au 
pire, la matrice de covariance est une matrice diagonale avec des termes 

35 infinis ou tres grands (de fagon a avoir des ecarts-type tres grands devant le 
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domains dans iequef le vecteur d'etat est susceptible d'evoluer) et Pestimee 
du vecteur d'etat le vecteur nu I, iorsque Ton n'a aucune idee sur les valeurs 
initiates. 

En pratique, le gain de correction d'un filtre de Kalman est 

5 "proportionnel" a Pincertitude sur Pestimation a priori de Petat du systeme et 
"inversement proportion ne!" a Pincertitude sur ia mesure. Si la mesure est 
tres incertaine et ['estimation de Petat du systeme relativement precise, Pecart 
constate est principalement du au bruit de mesure et la correction qui en 
decoule doit etre faible. Au contraire, si Pincertitude sur la mesure est faible 

10 et celle sur ('estimation de Petat du systeme grande, Pecart constate est tres 
representatif de Perreur d'estimation et doit conduire a une correction forte. 
C'est le comportement vers lequel on tend avec le filtre de Kalman. 

En resume, un filtre de Kalman sans vecteur de controle est defini 
par une matrice Fk correspondant a Pequation devolution definissant 

15 Pevoiution du vecteur d'etat du systeme physique modelise, une matrice H K 
correspondant a Pequation d'observation definissant les relations permettant 
de passer du vecteur d'etat au vecteur de mesure et une matrice de gain Kk 
mise a jour a Paide d'un processus iteratif mettant en jeu la matrice de 
covariance du vecteur d'etat P elle-mSme mise a jour au cours du processus 

20 iteratif et des matrices de covariance Qk et R K des bruits d'etat et de mesure. 

Le filtre de Kalman 3 opere au niveau des ecarts constates entre 
deux versions des informations de meme type provenant Pune, de la centrale 
inertielle INS 1 et Pautre, du recepteur GNSS 2, car cela permet de travailler 
sur des variables ayant des domaines de variation plus restraints sur 

25 lesquels ('approximation lineaire peut etre utilisee pour simplifier les 
equations de moderation et devolution (on parle alors de filtre aux erreurs), 
Dans ce cadre, Pequation devolution du filtre de Kalman modelise 
Pevoiution previsible des erreurs sur le systeme hybride lies aux residus des 
defauts capteurs inertiels non estimes et a la mecanisation de la plate forme 

30 inertielle (ici boucle fermee) de la centrale inertielle 1 qu'il deduit des 
differences constatees entre la prevision d'ecart qu'il fait et les ecarts 
effectivement mesures entre la position et ia Vitesse issues de la centrale 
inertielle et du recepteur GNSS. Elle est determinee par la matrice F dont la 
definition est fonction du mouvement du porteur, c'est-a-dire des parametres 

35 de vol dans le cas ou le porteur equipe du systeme hybride de 
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positionnement est un aeronef. La definition des diffe rentes versions de cette 
matrice F en fonction des parametres du mouvement du porteur sort du 
domaine de la presente invention. Elle ne sera pas detaillee dans la suite car 
elle est bien connue de I'homme du metier specialise dans le domaine de 
5 l'inertie. 

Les residus x K/K des erreurs du systeme hybride estimes a 
posteriori par le filtre de Kalman d'hybridation 3 sont utilises pour le recalage 
de la centrale inertielle INS 1 (boucle fermee). Les informations inertielles 
INS 1 etant regulierement recaiees, les points de position et de vitesse 
10 hybrides sont supposes plus fiables a condition que les mesures GNSS 
soient integres. 

Pour resumer, le traitement effectue par le filtre de Kalman 
d'hybridation se divise en trois temps : 

Propagation du vecteur d'etat et de la matrice de variance-covariance 
15 associee 

Xk*uk = F K+1 .X K/K + COR_FK__p K 

Recaiage du vecteur d'etat et de la matrice de variance-covariance a 
I'aide du gain Kk+i 

20 K K+1 = Pk+uk-H K+f {{{ K ^ .P K+U K H K+\ +R K+l) 

Pk+i/k+i - Pk+\/k ~K K+l .H K ^.P K+UK 
Application des corrections issues du filtre de Kalman 

d'hybridation a la centrale inertielle 

25 COR _FK_ p K = -X KlK 

La matrice d'observation H K +i permet I'observation des differences 
entre les positions et vitesses issues de la plate-forme virtuelle PFV et les 
positions et vitesses GNSS resolues. La matrice devolution F K+ i est calculee 
en utilisant les donnees de la plate-forme virtuelle PFV corrigees par le filtre 

30 de Kalman et permet la propagation de I'etat et de la matrice de variance- 
covariance associee. 

Le probleme rencontre avec les' systemes hybrides de 
positionnement INS/GNSS en boucle fermee est le risque de voir les 
operations de recalage propager a la centrale inertielle 1 les erreurs de 

35 position et de vitesse commises par le recepteur GNSS 2 en raison 
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d'informations defectueuses lui parvenant d'un ou plusieurs satellites 
defaillants (mauvaises corrections appliquees par le fiitre de Kaiman a la 
plate-forme virtuelle PFV, mauvaise estimation des defauts inertiels). 
L'integrite du systeme hybride de positionnement INS/GNSS en boucle 
fermee est alors difficile a assurer. 

La figure 3 illustre le schema d'un systeme hybride INS/GNSS en 
boucle fermee equipe d'un detecteur de probleme satellite operant a partir 
des points de position, en coordonnees geographiques, effeetues par le 
recepteur GNSS 2 avec tous les satellites visibles supposes au nombre de N 
et avec tous les satellites visibles moins fun d'entre eux, le satellite visible 
ecarte etant I'un quelconque d'entre eux. 

Le detecteur de probleme satellite 4 comporte un banc de N filtres 
predicteurs-estimateurs 40 1( 40 2 ,..., 40,,...,40 N associe a un banc de circuits 
de test 41 1f 41 2l ..., 41j,...,41 N controlant, par un circuit inhibiteur 5 intercale a 
la sortie du fiitre de Kaiman d'hybridation 3, la possibility d'un recalage de la 
centrale inertielle 1 . 

Les filtres predicteurs-estimateurs 40-,, 40 2 ,..., 40 i) ...,40 N 
fonctionnent en boucle ouverte. 

Le vecteur d'etat du ieme fiitre predicteur-estimateur 40j est 
compose des erreurs par rapport au vecteur d'etat du fiitre de Kaiman 
d'hybridation induites par la non prise en compte dans le recepteur GNSS 2 
des informations donnees par le ieme satellite visible et meme plus 
generalement par rapport a des mesures faites par d'autres equipements de 
navigation comme un altimetre barometrique ou un anemometre. 



erreurs _ attitude _ par _ rapport _ FK 
erreurs _ vitesse _ par _ rapport _ FK 
erreurs _ position _ par _ rapportFK 
erreurs _ biais _ aecelero _ par _ rapport _ FK 
erreurs _ derives _ gyro _ par _ rapprot __ FK 
erreurs _ baro _ par _ rapport _ FK 
erreurs _ anemo par _ rapport FK 



Err 1 - 
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Les matrices devolution F K des differents filtres predicteurs- 
estimateurs 40i, 40 2 ,..., 40j,...,40 N sont identiques a celle F« du filtre de 
Kalman d'hybridation 3 : 



Err l K + VK = F K+v Err^ /K 



(1) 



5 La mesure Z 1 d'un ieme filtre predicteur-estimateur 40 s est 

constitute de I'ecart entre la position, en coordonnees geographiques 
(latitude, longitude et altitude), resultant du point de positionnement P Nlot 
effectue par le recepteur GNSS 2 en tenant compte de tous les N satellites 
visibles et la position, en coordonnees geographiques, resultant du point de 

o positionnement P {N _ l)n effectue par le recepteur GNSS 2 en ecartant les 
informations envoyees par le ieme satellite : 



Z l = 



Ay' 

Az'' 



Nlot 



Les matrices H d'observation des differents filtres predicteurs- 

estimateurs 40-t, 40 2 ,.-., 40j 40 N sont toutes identiques et correspondent a 

I'equation d'observation permettant de passer du vecteur d'etat Err' au 

vecteur de mesure Z' . 

Les gains des differents filtres predicteurs-estimateurs 40i, 40 2 ,..., 
20 40j,...,40 N sont tous pris egaux a celui K du filtre de Kalman d'hybridation 3 
de fagon que le recalage sur les mesures z' qu'ils regoivent se passe comme 
le recalage du filtre de Kalman d'hybridation 3. 

Compte tenu de ces choix, les recalages des vecteurs d'etat Err' 
des filtres predicteurs-estimateurs 40 1f 40 2 ,.--, 40 it ...,40 N repondent a la 

25 relation : . 

Er4 +UK+1 = Er4 +1/K + K K+l \Z K J - H K+J .Err l K+UK ) (2) 

Pour detecter la defaillance de I'un au moins des N satellites 
visibles par le recepteur de positionnement GNSS 2, on effectue des tests 
statistiques sur les vecteurs d'etat £Vr'disponibles en sortie des differents 
30 filtres predicteurs-estimateurs 40i, 40 2 ,..., 40j,...,40 N au moyen du banc des 

circuits de test 41i, 41 2 ,..., 41 41 N . Ces circuit de test realisent des tests 

statistiques consistant a admettre une panne de satellite si, pour i allant de 1 
au nombre N de satellites visibles, Tune des deux inegalites sur les erreurs 
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de position des filtres predicteurs-estimateurs 40-,, 40 2 ,..., 40 ,40 N 

suivantes est verifiee : 



La variable K_seuil controle indirectement la valeur de I'erreur 
radiale acceptable. Elle est choisie en fonction du taux recherche de fausse 
10 alarme. Pour une probability de fausse alarme de 10" 6 , on prendra une valeur 
de 5,06 environ (repartition gaussienne) pour la variable K_seuil. 

La realisation de ces tests statistiques implique que I'on soit 
capable de propager et de recaler les variances associees aux vecteurs 
d'etat ErH des differents filtres predicteurs-estimateurs 40-,, 40 2 ,..., 40j,...,40 N . 
15 La propagation et le recalage peuvent se derouler ainsi : 




_ seuil.CO V _ Err^ tK (err _ lat) 



5 



ou 



\ Err Lv K (err _ Ion)" .Err l K+VK (err _ lon)\ > K seuiV X _ seuil.CO V 



ErrLuK&r Jon) 



A Initialisation , on pose : 



20 



COV _Err' /0 =[0] 

[o] etant la matrice nulle. 

A la premiere propagation, par definition : 




ce qui s'ecrit en raison de la relation (1) : 



25 




u 1 etant un bruit blanc (bruit d'etat). 
D'ou : 



30 



COV_E^ /0 =F l .COV_Err^ 0 .F l T 



F etant la matrice devolution du filtre de Kalman d'hybridation 3 et la 
matrice de covariance des bruits d'etat du filtre de Kalman d'hybridation 3. 
Au premier recalage, par definition : 
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15 



COV_Err' n = E^rri n JSrr^ T \ 
ce qui s'ecrit, en tenant compte de la relation (2) : 

En developpant cette expression, on trouve : 

...-4^/o .Errlj]H* .K X T -K v H v E[ErrU .Err^j] 



D'ou : 

COr_.Br/,, =(l-K 1 .H l ).COV_Err; i0 .(l-K ] .Hj + K l .E[z i Z iT \k 



1 etant la matrice identite. 

Le premier terme de I'expression du recalage de la matrice de 
covariance se calcule a partir de la matrice d'observation H 1 des filtres 

20 predicteurs-estimateurs 40 1s 40 2) ..., 40i,...,40 N et de la matrice de gain K, du 
filtre de Kalman d'hybridation 3. Le deuxieme terme se calcule au moyen des 
matrices de variance-covariance utilisees dans les calculs des moindres 
carres effectues pour la resolution des position, Vitesse et temps dans le 
recepteur GNSS 2, lors des points de positionnement a N et (N-1) satellites. 

25 En effet, lors de cette resolution, la methode des moindres carres evalue le 
degre de confiance sur les points resolus par r intermediate des matrices de 
variances covariances (matrices E[P (n -i)/i.P(n-i)/] pour les N points a (N-1) 
satellites et E[P lo t.Ptot T ] pour le point a N satellites). 

Or : , 



30 

Car : 
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20 



z 1 = 



Ay' 
Az' 



= P< 



En developpant la premiere expression, on obtient : 

Les deux premiers termes sont connus. lis sont directement 
fournis par les calculs des moindres carres lors de la resolution des differents 
points de position, Vitesse et temps effectues par le recepteur GNSS 2 

Les termes E[p [N _ l)n ^ t ] e t ^Jg^J, traduisant Hnter- 
correlation entre le point de positionnement P tot a N satellites et les points de 
positionnement P (N _ m a (N-1) satellites, peuvent etre evalues en reprenant 
I'expression des moindres carres de chaque point. 

En effet, on a : 

h In-i)h etant la matrice pseudo-inverse utilisee dans la resolution position, 
Vitesse et temps faite lors du point de positionnement effectue a N-1 
satellites (absence du ieme satellite) : 

h (n-i)/< =( H {N-i)/>- H (N-i)t) '.^^etant construite a partir de la matrice 
H {N-x)i<W\ est la matrice des cosinus directeurs liant les pseudo-distances d a 
la position resolue par la relation : d =H {lf _ 1)lt .P {N _ x)n 

d etant le vecteur forme par les pseudo-distances. 
Ainsi : 



25 Or : 



et 



E [P{ N -i)nPLh E 
= E 



d-d T ]H{ N _ l)n .HZ 



{ni-H lot y) T =(Hi.H tot y 



30 car la relation 
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{H T tor H tot y{HlH tot )=I 

entratne : 

<=> Hf Q{ .H m . -H toi ) } ~ / 

De ce fait : 

5 H{N-\)ii-H*l -iH{N-i)if^{N-i)ii) -H (at-i)/ i H tot-H tot) 

En outre, on a I'egalite suivante : 

10 La matrice H (N _ l)n ne differe de la matrice H tot que par la ieme 

ligne qui est nulle dans le cas de la matrice H {N _ 1)n et qui comporte les 
cosinus directeurs de la ligne a vue du ieme satellite dans le cas de la 
matrice H tot . En effectuant le produit de matrices Hl t .H {N _ l)n , les termes de 
la ieme colonne de la matrice Hj ol sont multiplies par les termes nuls de la 

15 matrice H {N _ )n . Ainsi, si les termes de la ieme colonne d'une matrice 
multipliee par la matrice H {N _ )n sont deja nuls (comme dans le cas de la 
matrice H( N _ l)n ), il n'y a pas de changement. Cela explique I'egalite ci- 
dessus. Par consequent : 

20 ={HL- H tot) 



25 



IT* r r *T 



done : 



Finalement, nous avons la relation : 

E[r.z iT \=E{p {N _ l)n .P^ )l ]-E\p ur pL\ 

qui montre que la matrice e[z 1 .Z' t \ peut etre evaluee a partir des matrices 
issues des calculs des differents points de position, vitesse et temps 
effectues par le recepteur GNSS 2. 
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Le premier recalage et la premiere propagation ne posent done 
pas de probleme. Les matrices devolution F-, et de gain utilisees sont 
celles du filtre de Kalman d'hybridation 3 tandis que la matrice d'observation 
H est commune a tous les filtres predicteurs-estimateurs 40-,, 40 2 
5 40i,...,40 N . 

Par recurrence, il est possible de montrer qu'il en est de me me 
pour les recalages et propagations suivants. 

Des qu'au moins un des tests statistiques est leve ou positif, il y a 
suspicion de mauvais fonctionnement de I'un des satellites visibles sur 
10 lesquels se base le recepteur GNSS 2 pour determiner position et vitesse du 
porteur. Ce franchissement est utilise pour declencher le circuit inhibiteur 5 
intercale entre la sortie du filtre de Kalman d'hybridation 3 et I'entree de 
recalage de la centrale inertielle 1, le satellite fournissant des donnees 
erronees ne sera reintegre que lorsque tous les tests statistiques ne seront 
15 plus leves, avec eventuellement un delai supplemental de precaution. 

On remarque qu'il est egalement possible d'estimer un rayon de 
protection sur la position hybride issue de la centrale inertielle 1 en utilisant la 
methode du maximum de separation. 

Le circuit 4 de detection de probleme satellite ne permet pas 
20 d'identifier le satellite fautif car ('information erronee provenant de ce satellite 
est prise en compte dans le point a N satellites visibles qui sert de reference. 

Sachant que I'un des satellites visibles utilises par le recepteur 
GNSS 2 ernet des informations suspectes, il est interessant de pouvoir 
('identifier pour le faire mettre a I'ecart par le recepteur GNSS 2. 
25 La f 'gure 4 donne un exemple de schema d'un systeme hybride 

INS/GNSS en boucle fermee equipe d'un detecteur de probleme satellite et 
d'un identificateur de satellite defaillant. 

Comme dans le cas de la figure 3, on retrouve une centrale 
inertielle INS 1, un recepteur GNSS 2, un filtre de Kalman d'hybridation 3 
30 operant sur les residus des erreurs entre la position, en coordonnees 
geographiques, donnee par la centrale inertielle INS 1 et celle, egalement en 
coordonnees geographiques, donnee par le recepteur GNSS 2 et fournissant 
des informations de recalage a la centrale inertielle INS 1, et un circuit 4 de 
detection de probleme de satellites analogue a celui de la figure 3 
35 commandant un circuit d'inhibition 5 intercale entre la sortie du filtre de 
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Kalman d'hybridation 3 et une entree de recalage de la centrale inertielle INS 
1. 

En plus de ces elements apparaTt un circuit 6 identificateur de 
satellite defaillant operant a partir de differents points de position, en 
5 coordonnees geographiques, effectues par le recepteur GNSS 2 avec tous 
les satellites visibles moins un et avec tous les satellites visibles moins deux, 
les deux satellites visibles ecartes etant a chaque fois differents. Ce circuit 6 
identificateur de satellite defaillant fonctionne en parallele du circuite 4 de 
detection de probleme satellite. I! est active par Ie circuit 4 en cas de 
10 detection de probleme de satellite. II identifie le satellite defaillant pour le 
recepteur GNSS 2 qui ne le prend alors plus en compte pour faire le point de 
position et de vitesse qui est envoye au filtre de Kalman d'hybridation. Une 
fois le satellite defaillant ecarte, le circuit 6 identificateur de satellite desinhibe 
le circuit d'inhibition 5 et la plate-forme virtuelie PDV de la centrale inertielle 
15 INS 1 est a nouveau corrigee. 

Le circuit 6 identificateur de satellite defaillant se base sur des 
tests en prenant comme references les differents points de position possibles 
avec N-1 satellites qui sont compares a des points de position a N-2 
satellites issus de ces points a N-1 . Ces tests consistent a : 
20 - evatuer, pour chaque selection possible de N-1 satellites parmi 

les N observables, les ecarts de position, en coordonnees 
geographiques, constates entre le point de position fait par le 
recepteur GNSS 2 avec I'ensemble des N-1 satellites de la 
selection et les differents points de position possibles a N-2 
25 satellites que peut faire le recepteur GNSS 2 en ecartant un 

satellite supplemental de la selection, 

soumettre ces ecarts a un banc de filtres predicteurs- 
estimateurs pour evaluer l'erreur induite conjointement par deux 
satellites, 

30 - comparer les etats de ce deuxieme banc de filtres predicteurs- 

estimateurs avec les variances qui leurs sont associees, 
identifier la famille de filtres ne levant pas leur test statistique et 
en deduire le satellite en panne, 

indiquer au recepteur GNSS le satellite corrompu pour le retirer 
35 de I'ensemble des satellites visibles utilises pour resoudre le 
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point de position (et de vitesse) a N satellites communique au 
filtre de Kalman 

desinhiber les corrections fournies par le filtre de Kalman a la 
centrale a inertie une fois le satellite defecteux retire du point du 
recepteur GNSS et le filtre reconfigure. 
En effet, ce n'est que lorsque le satellite fautif est evince de la 
selection de N-1 satellites et des selections de N-2 satellites qui en decoulent 
que les tests statistiques sur les etats des filtres predicteurs-estimateurs ne 
se declencheront pas (tests non leves). 

Comme montre a la figure 5, le circuit 6 identificateur de satellite 
defaillant comporte un banc de N.(N-1)/2 filtres predicteurs-estimateurs 61 u 
(1 <£<j:=M, i et j designent les satellites non observes lors de la resolution du 
point a N-2 fourni au filtre predicteur-estimateur) traitant toutes les 
possibility d'ecarts entre un point de position, en coordonnees 
geographiques, a N-1 satellites visibles et un point de position, egalement en 
coordonnees geographiques, a N-2 satellites visibles, le satellite exclu de la 
selection a N-1 satellite I'etant egalement de la selection a N-2 satellites. Un 
circuit de test 62 itj associe a chaque filtre predicteur-estimateur 61y du banc 
teste les etats modelises dans le filtre predicteur-estimateur par rapport a 
leurs variances. Un circuit de traitement 63 extrait des tests fournis par les 
circuits de tests 62 I(j I'identite du satellite fautif. 

Chacun des filtres predicteurs-estimateurs 61 (j du banc de filtre 
dedte a ('identification du satellite en panne a la meme configuration que les 
filtres predicteurs-estimateurs 40j du banc de filtres dedies a la detection de 
defaillance d'un satellite avec les rnemes matrices de gain K, devolution F et 
d'observation H, et recoit, comme mesure, I'ecart fourni par un circuit 
soustracteur 60ij entre la position, en coordonnees geographiques (latitude, 
longitude et altitude), resultant d'un point de positionnement P (AM)/ . effectue 
par le recepteur de positionnement par satellites GNSS 2 en tenant compte 
des N satellites visibles moins un, le i 6me et la position en coordonnees 
geographiques resultant d'un point de positionnement P {N _ 2)/U effectue par le 
recepteur de positionnement par satellites GNSS 2 en ecartant les 
informations du satellite deja exclu, le i §me et d'un autre, lef me . 

Les calculs des variances associees necessaires aux circuits de 
tests 62ij pour comparer les etats par rapport a leurs variances sont 
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exactement les mernes que ceux decrits pour les circuits de test 41 j associes 
aux filtres predicteurs-estimateurs 40j du premier banc et ne seront pas 
repris ici. Seule I'observation est differente puisqu'il s'agit ici de difference 
entre des points resolus avec (N-1) satellites visibles et des points resolus 

5 avec (N-2) satellites visibles. 

Le circuit de traitement 63 compare les tests des N differentes 
families, ou le f me satellite a ete enleve pour la resolution du point de position 
a N-1 et des points de position a N-2 satellites (1 ^^J), puis identifie le 
satellite defectueux par le fait que seuis les filtres predicteurs-estimateurs 

10 d'une familie, celle ne prenant jamais en compte le satellite defectueux n'a 
aucun test statistique leve. En effet, outes les autres families de filtres 
predicteurs-estimateurs comparant les points de position a N-1ou le satellite 
defaillant est present avec des points de position a N-2 satellites ont au 
moins un test statistique leve. Ce circuit de traitement 63 permet d'eviter 

15 ('utilisation du satellite corrompu dans la resolution du point GNSS foumi au 

systeme hybride. 

Pour faciliter la comprehension, les fonctions realisees dans un 
dispositif de surveillance de Tintegrite d'un systeme hybride INS/GNSS ont 
ete illustrees sous forme de blocs separes mais il est bien evident qu'elles 
20 peuvent etre remplies par un meme calculateur a logique programme, par 
exemple le calculateur de gestion du vol si le porteur du systeme hybride de 
positionnement est un aeronef. 
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REVENDICATIONS 

LDispositif de surveillance de I'integrite d'un systeme hybride 
constitue d'une centrale inertielle INS (1), d'un recepteur de position nement 
par satellites GNSS (2) operant a partir d'une constellation de N satellites 
5 visibles, et d'un filtre de Kalman d'hybridation (3) ayant un vecteur d'etats 
correspondant aux erreurs du systeme hybride, en particulier les erreurs de 
biais et de derive residuelles de la centrale inertielle INS (1), observant les 
ecarts entre les points de positionnement et de vitesse fournis, en 
coordonnees geographiques, par la centrale inertielle INS (1) et par le 

10 recepteur GNSS (2), ayant une matrice devolution F modelisant revolution 
des erreurs du systeme hybride, une matrice d'observation H modelisant les 
relations entre le vecteur d'etats et les ecarts observes entre les positions et 
de vitesses delivrees par la centrale inertielle INS (1) et le recepteur GNSS 
(2), et un gain K minimisant la covariance de I'erreur faite sur I'estimation a 

15 posteriori du vecteur d'etats du filtre de Kalman et en particulier des erreurs 
residuelles de la centrale inertielle INS (1), et delivrant une estimee a 
posteriori des erreurs du systeme hybride qui sert a recaler la centrale 
inertielle INS (1), caracterise en ce que le recepteur GNSS (2) delivre en plus 
d'un point de position effectue a partir des N satellites visibles qu'il a en vue, 

20 des points de position P (N . 1)/i , i e [1...N] resolus avec N-1 satellites visibles 
deduits de la constellation des N satellites visibles en la privant a chaque fois 
d'un satellite different et en ce qu'il comporte un circuit (4) detecteur de 
probleme satellite comportant un banc de N filtres predicteurs-estimateurs 
(40i) de I'erreur induite par le satellite qui a ete 6te lors de la resolution d'un 

25 point de position a N-1 satellites, ayant le gain K et la matrice devolution F 
du filtre de Kalman d'hybridation (3), observant chacun I'ecart entre le point 
de position, en coordonnees geographiques, delivre par le recepteur GNSS 
(2) en observant les N satellites visibles et I'un des points de position P ( N-i)/i, 
egalement en coordonnees geographiques, delivres par le recepteur GNSS 

30 (2) en observant N-1 satellites visibles et des circuits de test (41 0 comparant 
les etats des filtres predicteurs-estimateurs (40j) a leurs variances et 
detectant une panne satellite lorsque le test est positif, I'ecart trouve etant 
superieur a un seuil de detection. 
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2. Dispositif selon la revertdication 1, caracterise en ce que les 
seuils de detection sont des seuiis statistiques qui tiennent compte de la 
covariance associee au type d'erreur de positionnement considere. 

5 3. Dispositif selon la revendication 1, caracterise en ce que les 

seuils de detection sont des seuils statistiques qui tiennent compte de la 
covariance associee au type d'erreur de positionnement considere et qui 
sont fonction du taux de fausse alarme acceptable pour le test. 

10 4. Dispositif selon la revendication 1, caracterise en ce qu'il 

comporte un circuit (5) d'inhibition de recalage intercale entre la sortie du 
filtre de Kalman d'hybridation (3) et une entree de recalage de la centrale 
inertielle INS (1), et active par le circuit (4) detecteur de probleme satellite. 

15 5. Dispositif selon la revendication 1, caracterise en ce que les 

filtres predicteurs-estimateurs (40j) du circuit (4) detecteur de probleme 
satellite ont tous la meme matrice d'observation. 

6. Dispositif selon la revendication 1, caracterise en ce que le: 

20 recepteur GNSS (2) delivre des points de position P (N -i>/i i e [1 f ..N] et P ( N-2)/ij S 
(1 <l<j<N), en coordonnees geographiques, resolus a partir de N~1 et N-2 # 
satellites deduits des N satellites visibles en retirant a chaque fois un satellite # 
visible different, le f me pour le point P(n-d/i, et en retirant en plus un autre 
satellite parmi les N-1 restant, le f me pour le point P ( N-2)/i,j, et en ce qu T il 

25 comporte un circuit (6) identificateur de satellite defaillant comportant un 
banc de Nx(N-1)/2 filtres predicteurs-estimateurs (61 if j) estimateur de Terreur 
induite conjointement par deux satellites parmi les N observables ayant le 
gain K et la matrice devolution F du filtre de kalman d'hybridation (3), 
observant chacun Tecart entre un point de position (P(N-i)/0, en coordonnees 

30 geographiques, delivre par le recepteur GNSS (2) a partir d'une constellation 
specif ique de N-1 satellites visibles et Tun des points de position P(N-2)/ij 
delivres par le recepteur GNSS (2) a partir d'une constellation de (N-2) 
satellites visibles deduite de la constellation specifique de (N-1) satellites 
visibles privee de Tun de ses satellites, le f me , les filtres predicteurs- 

35 estimateurs pouvant etre regroupes par families de N-2 elements en fonction 



1er depot 
26 



de la constellation specifique de N-1 satellites prise en compte, un banc de 
circuits de test (62 u ) comparant les etats des filtres predicteurs-estimateurs 
(61i,j) a leurs variances et detectant une anomalie lorsque I'ecart trouve est 
superieur a un seuil de detection, et un circuit de traitement des resultats des 
tests (63) centralisant les tests leves identifiant, en cas d'anomalie detectee, 
un satellite fautif comme etant le satellite exclu d'un point de position resolu' 
avec N-1 satellites dont la famille de filtres predicteurs-estimateurs est la 
seule a n'avoir aucun de ses elements detectant une anomalie. 

7. Dispositif selon la revendication 6, caracterise en ce que les 
seuils de detection des circuits de test (62 u ) du circuit (6) identificateur de 
satellite defaillant sont des seuils statistiques tenant compte de la covariance 
associee au type d'erreur de positionnement considere. 

8. Dispositif selon la revendication 7, caracterise en ce que les 
seuils de detection des circuits de test (62 rj ) du circuit (6) identificateur de 
satellite defaillant sont des seuils statistiques qui tiennent compte de la 
covariance associee au type d'erreur de positionnement considere et qui 
sont fonction du taux de fausse alarme acceptable pour le test. 

9. Dispositif selon la revendication 6, caracterise en ce que les 
filtres predicteurs-estimateurs (61 y) du circuit (6) identificateur de satellite 
defaillant ont tous la meme matrice d'observation. 

25 10 - Dispositif selon la revendication 6, caracterise en ce que les 

filtre predicteurs-estimateurs (40 ia 61 y) du circuit (4) detecteur de probleme 
satellite et du circuit (6) identificateur de satellites defaillant ont tous la meme 
matrice d'observation. 
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